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Abstract 

We study variants of one-dimensional g-color nearest-neighbor voter 
models in discrete time. In addition to the usual voter model transitions 
in which a color is chosen from the left or right neighbor of a site there 
are two types of noisy transitions. One is bulk nucleation where a new 
random color is chosen. The other is boundary nucleation where a ran- 
dom color is chosen only if the two neighbors have distinct colors. We 
prove under a variety of conditions on q and the magnitudes of the two 
noise parameters that the system is ergodic, i.e., there is convergence to 
a unique invariant distribution. The methods are percolation-based us- 
ing the graphical representation of the model which consists of coalescing 
random walks combined with branching (boundary nucleation) and dying 
(bulk nucleation). 

1 Introduction 

In this paper we consider a class of one-dimensional interacting particle systems 
with a focus on ergodicity, i.e., whether there is always convergence to a unique 
invariant distribution. The evolving state of the system is an assignment of 
colors from {1, q} to the sites in Z and the transition probabilities of a given 
site depend only on the color of its left and right neighbors. We actually study 
a two-parameter, e and 5, family of models which includes classic voter models 
with or without noise and stochastic Ising and Potts models at zero and nonzero 
temperature. Together the two parameters control the rates of bulk nucleation 
of random colors and boundary nucleation when the two neighbors disagree in 
color. Our main results are about ergodicity under various conditions on e, 5, 
and q. The methods we use are percolation-based and rely on the fact that the 
dual system is a model of coalescing random walks that also have branching and 
dying. 

We start by recalling the classic discrete time voter model in one dimension 
|I3| . The voters or particles are located at each site of Z and take one of the 
q possible opinions (or colors). At each point in time t the particles make 
decisions on whether to keep or change their color. The decision consists of 



'Research supported in part by NSF grants OISE-0730136 and MPS-1007524 



1 



randomly choosing one of the two neighbors and adopting its color from time 
t — 1. The process of choosing a neighbor may be represented by drawing a 
random arrow to one of the nearest neighbors at time t—1 and thus the history 
or genealogy of the color of a particular site may be traced to time zero by the 
path of backward arrows coming out of each space-time site. The distribution 
of the genealogy paths is that of coalescing simple symmetric random walks. 

Noise can be introduced in the voter model by having each site make a 
decision as before with probability p or else choose a color uniformly at random 
out of the q possible ones with probability I — p. The genealogy of the noisy 
voter model is represented by the coalescing random walks with "dying" , where a 
path dies at a point where a color choice was made uniformly at random. Such 
points, which we call bulk nucleation points, represent death of a genealogy 
path or birth of an opinion. If the opinions are +1 or -1, the noisy voter model 
coincides with the stochastic Ising model with Glauber dynamics. 

The models we consider in this paper allow for the number of colors of a 
single site to be any positive integer q > 2. The decision that each particle 
makes at each point in time has the two possibilities for a noisy voter model: 
one is to choose a color of a randomly chosen neighbor and another is to gener- 
ate a uniformly random color. A third possibility is to take the color of the two 
neighbors if their colors agree or to generate a uniformly random color if they 
disagree (boundary nucleation). A site where such a decision is made we rep- 
resent by a double-arrow in the genealogy graph. To determine the probability 
of these three choices we introduce two parameters, S and e. This model, for 
certain choices of 6 and e, is related to stochastic Potts models with q colors (for 
the relation with zero temperature continuous time Potts models see [5], with 
continuous time nonzero temperature Ising models see [9] , with nonzero temper- 
ature Potts models in continuous and discrete time see [Hj). In the space-time 
coordinates the evolution of the system may be represented by coloring vertices 
of Z^. Since a color of vertex {z,t) depends on colors of vertices {z — l,t — 1) 
and (z + 1, t — 1) only, there are actually two independent systems evolving: one 
lives on the sub-lattice with z + t even and another with z + t odd. Thus, we 
can restrict our lattice to be the part of which contains only vertices with 
z + t even. 

The special case where the number of colors q -> oo is also considered. In 
this case it is convenient to view coloring as partitioning sites into equivalence 
classes of the same color rather than coloring the sites, and at nucleation points 
the color (or equivalence class) that is supposed to be assigned uniformly at 
random is chosen to be completely new and different from all the previously 
existing ones. 

The main issues that we investigate can be described by the term memory 
loss: to what extent does the system in the long run depend on the initial condi- 
tions. The existence of an invariant distribution is known for all lattice systems 
in discrete and continuous time and is based on compactness of the state space. 
See [7] for discrete time systems and [TB] for continuous time systems. If the 
invariant distribution is not unique, then memory is not lost. Limiting behavior 
will depend on the initial distribution as there are at least as many possible 
limits as there are invariant distributions. Ergodicity means convergence to 
a unique possible limit for all initial conditions and is the strongest degree of 
memory loss. We prove ergodicity under certain conditions on the parameters 
S and e. For others the question remains open. For all nonzero parameters we 
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prove a weaker type of memory loss - color permutation invariance of any limit 
or sub-sequence limit distribution. This shows that if one begins with all sites 
of the same color, any sub-sequence limit distribution does not depend on the 
starting color. 

There are natural diffusive scaling limit analogues of the model considered 
in this paper where space-time is continuous rather than discrete and random 
walks become Brownian motions. Some of these, such as the Brownian web 
and net and related continuum voter models, have already been studied [3] and 
others are in preparation |14j . One reason we focus here on the discrete time 
model is because the only existing published result [15] about convergence to 
these limits starts from the discrete time lattice model. Another reason is to 
allow us to compare ergodicity and percolation properties of the model - see 
the Remark following the proof of Theorem 2 in Section 3 below. We expect 
that the discrete model results of this paper can be extended to the continuum 
setting, but do not explore that here. We note however that there are a number 
of papers [H |4l |5l [10] that treat related continuous time lattice models that also 
have a natural percolation substructure to their dual models. 

A formal definition of the model and presentation of main results is given 
in the next section. Section 3 contains all the proofs. The proofs are based on 
the underlying directed arrow percolation process which traces the genealogy 
of colors. Percolation methods such as dynamic renormalization (see [3]) or 
enhancement (see |T]) that we use in this paper are described in detail in [TT]. 
The enhancement argument we utilize (see Lemma 2) is a close adaptation of 
that of Aizenman and Grimmett (PP or [H]). We present the argument in detail 
here because the percolation model that we are considering does not fall into 
the class of models considered in jlj: our model is directed and there are more 
than two possible states of the arrows at each site. But it does fall into the class 
of models treated in [T^]. 

2 Main Results 

Let V be the sub-lattice of with the sum of space and time coordinates even. 
At each vertex v define a random variable X{v) that takes values from the set 
of arrow configurations { \ , , \/^ , } with probabilities P{X{v) =\) 
= P{Xiv) = 1(1 - S), P(X(v) = 8{l - e), and P{X{v) = 0) = ^e, 

where 5 and e are parameters with values between and 1. The arrows coming 
out of a vertex v reach to the nearby vertices in the row above to the left in case 
X(w) = "\", to the right in case X(i;)= V" of w, or both in case X(v)="\/>". 
The outcome "0" means there are no arrows (bulk nucleation) . This will define 
our X-arrow percolation model. (Note that the genealogy paths of the noisy 
voter model can be recovered by setting e — 1 and for the classic voter model 
by setting 5 — 0.) At the vertices of the same lattice place color-valued i.i.d. 
random variables Y{v) with each color drawn uniformly out of q possible colors 
for some finite q > 2. To define the color process Z{v) take the bottom half of 
the lattice with the top row being the row containing the origin (0,0). Direct 
time (row count) down, opposite of the direction of the X-arrows. Values for 
the Z(vys will be constructed using initial (time-zero) values together with the 
X(w)'s and y(u)'s, as follows. At each vertex of the top (initial) row assign 
Z{v) colors for all vertices v according to some initial distribution. For each 
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subsequent (lower) row following the initial row determine the colors of each 
vertex v according to a rule that depends on the values of X(v), Y{v) and the 
Z values of the two adjacent vertices in the row just above of v called vi for the 
vertex above left of v and Vr for the vertex above right of v. The rule is the 
following 



Denote by Zn the random variable that represents the sequence of colors Z{v) 
of all the vertices of row n. Zn is a discrete time Markov process. Here we 
investigate under what conditions it is ergodic. The existence of an invariant 
distribution for Z„ is a known fact that is explained in [7] among many other 
sources. By ergodic we mean that the invariant distribution for Zn is unique 
and that any initial distribution converges to it. 

The first question of interest is whether X-arrows percolate. To say it more 
precisely, let a path of X-arrows be a sequence of vertices vi,V2,V3 ... such that 
Ui+i lies above right or above left of Vi with X{vi) containing the arrow from 
Vi to Vi-i-i. A finite path is said to terminate at Vn when X(z;„)="0". The X- 
arrows are said to percolate if there exists with strictly positive probability an 
infinite path of arrows starting at the origin or, equivalently, if with probability 
one there exists an infinite path of arrows starting from somewhere. 

Theorem 1. V (5 > 0, 3 edS) > such that when e > edS), the X-arrows do 
not percolate and when e < ec{5), the X-arrows percolate. 

As a corollary of Theorem 1 it is shown that when the X-arrows do not 
percolate the Zn process is ergodic for any value of q. Although we have not 
shown ergodicity for all positive e, the following theorem extends the interval of 
e for which the Zn process is ergodic to some positive distance below ec((5) where 
the X-arrows do percolate. Thus ergodicity is not equivalent to non-percolation. 

Theorem 2. V (5 > and q>2,3 e'^S) < ec(^) such that for e > e^((5), the Zn 
process is ergodic. 

A weaker result is established for all non zero e: 

Theorem 3. V(5 > 0, e > 0, 2 < q < oo, and any initial configuration Zq, any 
.sub-sequence limit distribution is color permutation invariant. 

This result does not imply convergence to a limiting distribution nor does 
it imply uniqueness of the invariant distribution as it is still possible to have 
different fc-dimensional marginal distributions for k > 2. 

There are two special cases where the Zn process is ergodic for all positive 
S and e. In the case q — 2, the ergodicity of Zn is known. This follows from 
the fact that for g = 2, at "\/>" points when the colors of neighbors disagree 
choosing color 1 or 2 uniformly at random is equivalent to choosing the left or 
right neighbor uniformly at random. Therefore, the paths of arrows that are 
used to trace the color genealogy of each vertex are distributed as coalescing 
simple symmetric random walks with dying ( "0" points). Such paths are almost 
surely finite, and ergodicity will follow by the argument of Corollary 1 below. 



Z{v) 



Z{vi) 

Z{Vr) 

Z{vi) 
Yiv) 



if X{v) = \, 
zf X{v) = ^, 

if X{v) and Z{vi) = Z{vr), 

if X{v) and Z{vi) ^ Z{vr) or X{v) = 0. 
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For q ~ CO, instead of colored vertices the state space of Z„ may be considered 
to be partitions of Z into infinitely many possible equivalence classes that play 
the role of colors. In this special case the ergodicity is established by Theorem 
4. 

Theorem 4. When q = oo, WS > and e > 0, the Zn process is ergodic. 

The conjecture of ergodicity for 2 < q < oo for all (5 > and any e > 
remains open. 



First we clarify a possible source of confusion: the time (row count) for paths 
of arrows is in the direction of the arrows which is upward, but the time for the 
color evolution process Z„ is against the arrows or downward. This is due to 
the color evolution rules that use arrows pointing to the colored vertices from 
the past in determining colors for the vertices in the present-i.e.,the arrows 
track the color genealogy. Throughout the rest of the paper we will choose the 
direction of time depending on whether we consider colors or not. 

Proof of Theorem 1: Fix (5 > 0. For simplicity let x describe the horizontal 
coordinate of the lattice V and y describe the vertical coordinate of the lattice. 
Let 6(e) be the probability of having an infinite path from the origin. The fact 
that 6(e) decreases as e increases follows from a coupling argument similar to 
the one used in the standard percolation model as we now explain. For each 
vertex v of the lattice let U{v) be a collection of i.i.d. uniform[0,l] random 
variables. For each e > set 



For a fixed e, the marginal distribution of X^{v) is the same as the distribution 
of the arrow- valued random variable X{v) with the parameters 5 and e. For 
ei < e2, has all the arrows of X^^. Hence, 6(ei) > 6(e2). To show that 
Cc < 1 we observe that the number of vertices at height n reached from the 
origin by paths of arrows is dominated by the number of vertices at height n of 
the branching process with the same distribution of offspring as the number of 
arrows X{v) coming out of a vertex v. It follows that since of offspring) = 
1 • (1 — (5) + 2 • (5(1 — e) + • (5e < 1 for e > -i then the branching process eventually 
almost surely dies out, which in turn implies that 6(e) = for such e. 

It remains to prove that Cc > 0. The idea of the proof (see Lemma 1 below) 
is taken from an argument of Durrett [5] giving an upper bound for the critical 
probability of the oriented independent site percolation model. It is shown in 
[5] that for oriented independent site percolation on Z^, pc < 0.819. We use a 
dynamic renormalization technique [3] . The idea behind it is described in . 
Define a box B(k,n) to be Qk wide and n high, where k and n will be chosen 
later. Such a box will be a renormalized site in an oriented independent site 
percolation model. Let the lower, upper, and side boundaries of the box be rows 
or columns of vertices parallel to and just inside the edges of the box. Position 
the boxes in rows such that the boxes of the first row have their lower boundary 



3 Proofs 



X,{v) 







if U{v) e [0,^), 
if U{v) e i^A-S), 
if U{v) e [1-^,1-^e), 
if U{v) e [l-^e,l]. 
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on the t ^ 1 line (here t goes up) with a column of vertices between adjacent 
boxes' side boundaries. The second row of boxes will have its lower boundary 
right above the upper boundary of the first row but shifted by 3A: + 1 to the side 
so that the middle of each box in the second row is above a column of vertices 
left in between boxes in the first row. The third row will be positioned on top 
of the second row similarly, and so on. See Fig. [T] for a diagram. 



row 3 



row 2 



row 1 



Figure 1: Relative positions of boxes. 



Since the boxes do not overlap, the events described in terms of the vertices 
inside of a box are independent for different boxes. For each box we define the 
following events. Divide the lowest row of vertices in the box into six equal parts 
with the length of each part equal to k and place the x — coordinate in the 
middle of the row. For every vertex v in the interval / — [—2k, 2k] of the lowest 
row of vertices of the box let be the event that there exist paths from v to at 
least one of the vertices in the interval Fi = [—3k, —k] and to at least one of the 
vertices in the interval F2 = [k, 3k] in the horizontal row of vertices right above 
the box (which belong to the boxes in the next row of boxes) . Both paths have 
to be completely inside the box containing v. Note two things: i) even though 
the paths reach one level above the box containing v, their existence depends 
on the arrow states of vertices inside the box only; ii) the destination intervals 
Fi and F2 are in corresponding intervals / of the two boxes in the row of boxes 
above to the left and to the right of the box containing v. For a diagram of the 
described event see Fig. [2] 

Lemma 1. // min^jg/ P(A„) > pc, the critical value for oriented independent 
site percolation on , then the X -arrows percolate. 

Proof of Lemma 1: The lemma follows by comparing boxes to sites in 
the oriented site percolation model using the standard dynamic renormalization 
technique described in [3] and [TT] among other places. The idea behind the 
technique is to construct the infinite cluster of the origin sequentially. Consider 
the box containing the origin occupied if Aq holds. Then in a predetermined 
order that increases with height check the boxes that are upper-right and upper- 
left neighbors of already occupied boxes. A box is checked if one of its neighbors 
below- left or below-right or both are occupied. By an inductive argument, this 
implies that the interval / of the box being checked contains at least one vertex 
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-3k -2k 2k 3k 

Figure 2: Box is occupied given it is checked. 



V with a path of X-arrows from the origin reaching it. Choosing one such v 
we declare the box occupied if Ay holds. Under the assumption of Lemma 1, 
the cluster of occupied boxes will stochastically dominate the cluster of vertices 
reached from the origin in the oriented independent site percolation model. □ 
To apply Lemma 1 we will choose appropriate values for the height of a box 
n and its width 6k that make the probability of the event Ay arbitrarily close 
to one for all possible starting vertices v in the interval I = [—2k, 2k]. First, 
set the parameter e to 0, which sets the distribution of the outgoing arrows for 
each vertex to P(\) = P(/') = ^(1 - P(\/") = 6, and P(0) = 0. Pick 
an arbitrarily small C > 0. Now we will determine dimensions n and 6k of a 
box that will make P{A-2k) above 1 — 5- Let v be the vertex at the left end of 
the interval I of the box (with horizontal coordinate —2k relative to the middle 
of the box and vertical position 1 relative to the bottom of the box). Note 
that with e set to zero, at each time step we equally likely advance either to 
the right or to the left and sometimes both. Let S\ denote a path starting at 
this V which follows the arrows and makes a random choice at a site with both 
arrows. If we consider the rightmost path S2 from v, we advance to the right 
with probability P(/') + P(\/') = ^(1 - 5) + ^ = i(l + 5) and we advance to 
the left with probability P(\) = ^(1 — 5). By translating v to the space-time 
origin we obtain coupled simple random walks Si and starting at the origin 
with 5i symmetric and 52 asymmetric. A corresponding translate of the event 
Ay then contains the intersection event {Si hits the interval [—k, k] at time n 
without going below —k} p| {<S'2 hits the interval [3A;,5fc] at time n without 
going above 5k}. Choose k = The probabilities of the events defined in 
terms of the random walks can be easily approximated using the central limit 
theorem: 

and the reflection principle: 



P{Si{i) > -k for alio <i<n) = l- 2P{Si{n) < -fc) w 1 - 2P{Z < - 
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where Z is a standard normal random variable. Similarly for 5*2, whose incre- 
ments Xi have E{Xi) — 5 and Var{Xi) = 1 — (5^, we have 



P{S2{n) e [3/c,5fc])«P(Ze [- 



\fnb 



]) 



and 



P{S2{i) <5fcV0<i<n)>l - 2P{S2{n) > 5fc) « 1 - 2P(Z > 



)• 



g by taking n large enough 
Once we have 



All these probabilities can be made greater than 1 
so that the probability of their intersection is at least 1 — § 
P(A_2fc), and by symmetry P{A2k), above 1 — ^, since there is at least one 
arrow coming out of every vertex, it is easy to see from Fig. [3] that for any 
V e [-2fc,2fc], P{Ay) > P(A_2fen^2/c) > 1-C- (As noted by M. Damron, since 
also PifXiZ^-2k — -P(^-2fe n ^2fe) > 1 — C: one can replace our dynamic by 
a static renormalization argument.) 




-3k -2k 



2k 3k 



Figure 3: A^2k n A2k C Ay for -2fc <v<2k. 



Now, as we have the size of a box fixed by choosing n large enough, the 
probabilities of the events defined in terms of the arrows inside of a specific box 
are polynomials in e and, therefore, continuous in e. If the values for P{Ay) for 
all the vertices v in [—2k, 2k] exceed 1 — C for e = 0, then they exceed 1 — 2C for 
some small e > by continuity. This completes the proof of Theorem 1. □ 



One consequence of Theorem 1 is that, for fixed 5 and e > edd), the cluster 
of arrows coming out of vertices of some finite interval of the initial row is almost 
surely finite. All the paths from those finitely many vertices terminate at "0" 
points. For the color process Z„ this means that with the initial time long 
enough in the past the distribution of colors of any finite interval at time for 
a fixed S and e > edS) is determined by the distribution of a finite cluster of 
arrows coming out of the interval and by the i.i.d. Y{v) colors at the end ("0") 
points of the cluster. This proves the following corollary. 

Corollary 1. The Z„ process is ergodic for any 5 > Q and e > tc{5)- 
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To prove Theorem 2 we will use an auxiliary lemma which is an adaptation 
of the enhancement results of [11] for standard bond percolation. The lemma is 
for the arrow configurations and does not use any color-valued variables. First 
we describe the types of enhancement we will use. On the lattice V let C„ be 
the deterministic subset of V which consists of all the vertices up to height n 
that can be reached from the origin by arrows taken from the set of all possible 
arrow configurations {\, /^,\/^, } at each vertex v (see Fig. 0]). 



Figure 4: The sets of vertices CI, C2 and C3. 



Let n = {\, }^ , and A = {0, 1}^. For any v £ V and uj £ n, 

let F{uj{v)) be an enhancement function, where i-u(v) is a restriction of uj to the 
subgraph v + Ci. F{uj{v)) is defined as follows: if the arrow value at v is 
"\/"' and at least one of the two remaining vertices in v + Ci has value "0" , 
then the arrow value at v is changed to "0" while arrow values at other sites 
are unchanged; in any other case no change is made — see Fig. [5] 




Figure 5: Configurations changed by enhancement function. 



For any configuration uj £ Q, and A e A, at each vertex v £ V , the en- 
hancement Lo{v) -> F{ijj{v)) is made if X{v) = 1. Fix 5 > Q. The probability 
measure we take on D. is the product measure joint distribution of the previously 
described i.i.d. X{v) random variables at the vertices v £V . The probability 
measure we take on A is the product measure joint distribution of i.i.d. Bernoulli 
random variables at the vertices v with probability of "1" equal to s and prob- 
ability of "0" equal to 1 — s. The probability measure on $7 x A denoted by P^ s 
is the product of these two measures on VL and A. Let G'^"''^{uj, A) be the graph 
of arrows from the space 17 x A that is obtained from the random graph G{uj) 
of arrows of Q (whose marginal distribution is the same as the original arrow 
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process distribution) by applying the enhancement function F{ijj{v)) at each v 
with X{v) = 1. Define 

^e7ih _ j^fjg . p^,,{Q'^'^^(ijj^ A) has an infinite path from the origin) = 0} 
and 

tc = inf{e ; P^^s{G{uj) has an infinite path from the origin) — 0}. 
We are ready to state our auxiliary result: 
Lemma 2. For any s > 0, e^"'' < Cc- 

Proof of Lemma 2: In the space $7 x A define the event An={There is 
a path of arrows in G^^^{uj,X) from to the top of C„}. Let 0„(e, s) = 
Pn^si^n)- For any v € V, the event {v is uj — pivotal for An} is the collection 
of aU {lu, A) such that for w = w' = w" off v, "\/" and "0", 

Ia„{oj',X) ^ Ia^{oj",X) . Here, Ia is the indicator function of the event A. 
It should be noted from the definition above that the event {v is uj — pivotal 
for A„} does not depend on the value of lu at the vertex v itself. Similarly, 
the event {v is X — pivotal for An} is the collection of all (w. A) such that 
/a„(w. A') 7^ IaA^, A") with A = A' = A" off v, X'{v) = 0, and X"{v) = 1. 

In the remainder of the proof we show that there exists a path from (e', s) 
to (e",0) (with e' < ec and e" > ec) along which 6„(e, s) is nondecreasing for 
all n. First we adapt the proof of Russo's formula from [TT]. Set N = |C„| 
and denote the vertices in C„ by w^, 1 < i < Suppose that for every vertex 
Vi G Cn the value of e used to describe the distribution of X(vi) is a distinct 
variable e^. For each Vi, let Ui be a uniform[0,l] random variable with the Ui^s 
independent for distinct i. If we assign 



uj{vi) 



if U, e [0,^2^), 
if a, e [^,i~S), 
if U, e[i-5,i- 5e, 
i/ f/, e [l-(5e„l]. 



the new model will correspond to the model on 51 described above (when = e 
Vi). Then 

9e„(ei, EN, s) e„(ei, ei + h, e^^s) - e„(ei, e^, eat, s) 
^ = 1™ 1 

—ShPf: s(vi is UJ — pivotal for An) 
= lim ■ ; 

h^O h 

= —SP^ s(vi is UJ ~ pivotal for An). 

The second equality follows from the fact that configurations that are counted 
in 0„(ei, Ei, e^v, s) but not in G„(ei, + /i, ejv, s) are those where 
every path to the top of C„ from the origin must use one of the two arrows 
coming out of v, and the value of Ui is between 1 — S{ei + h) and 1 — Ssi. To 
compute the value of ^^'g^^^'^-* we use the chain rule and get 

^'^"g^' = -S^^P^^sivi is UJ - pivotal for An). (1) 
1=1 
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A similar argument shows that 



90 (e s) ^ 

^— ^ — = - ^ Pe,s{vi is A - pivotal for An). (2) 

i—l 

A vertex v is ^'uj — pivotal for An" if and only if one of the following two disjoint 
cases holds: 

i) every path from the origin to height n goes through w; 

ii) the only paths from the origin to height n that do not go through v go 
through a vertex u{v) below left or below right oiv such that (a) , 
(6) X{u{v))=l and (c) making the value of oj{v) to be "0" allows for the en- 
hancement to change the value of uj{u{v)) to "©" killing all the paths. 

I II III 

level n 



Figure 6: Examples of the configurations in the events I, II, and III. 



We next decompose case i) into two disjoint events that we will denote by 
/ and II; case ii) will be denoted ///. Event / is when the value of Ld{v) fixed 
to be "\/^" allows for the enhancement to change it to "0" . As seen in Fig|6l 
event / contains configurations where all possible paths from the origin to the 
initial level come to v and continue from one of its upper-left or upper-right 
neighbors having the other one at state "0" . Whether the paths exist depends, 
of course, on the state of uj{v), but whether v is w-pivotal does not (see the 
definition of w-pivotal in the beginning of this proof). However, for v to be uj- 
pivotal X(v) must be 0, or otherwise the possible paths would not exist whether 
uj{v) is "0" (obviously) or (would be killed by enhancement). Then by 

relaxing the requirement that A(w)="0" and by fixing uj(v) to be '"X/"', we 
make v also "A — pivotal for An' ■ 

Event // is when the enhancement leaves the value of aj(w) fixed to be 
unchanged (if the children of v both have outgoing edges); then by switching 
the value of ui at the vertex above left or above right that is not involved in 
connecting the origin to height n (if both are involved pick the left one) to "0" we 
make v to be '^X — pivotal for An" ■ Therefore 6{l — e)j^P^g(^I) < Pe,s{v is A — 
pivotal for An) and (5(1 — e)deP^^s{II) £ Pe.siv is X — pivotal for An). 

In event /// (i.e. in case ii)), by relaxing the requirement that X{u{v)) = 1 
(if both below right of v and below left of v vertices are such as described in 
case ii), then pick u{v) to be the left one) and by fixing uj{v) to be "0" we make 
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u{v) to be "A — pivotal for An". Therefore, ^6eP^^s{III) < Pe.s{u{v) is A 
pivotal for An). Now 

P,.s{v ISLU- pivotal for An) = PeAl) + Pe.sill) + PeAlU) 

1 — s 

< ^Pe,s(w is A - pivotal for An) 

+ Tn- — 7T"^>s,s(w is \- pivotal for An) 
o(l — ej de 

+ s-^[PtA'^beiow ieft{v) is A - pivotal for An) 

+ Pe,siubelow right{v) is X ~ pivotal for An)]. 

By summing over all the vertices in C„ we get 

JV 

^^Pe^i(t;i is Lo — pivotal for An) 

i=l 

1 1 ^ 

< ~ * + ^) 51 P>i.,s{vi is X- pivotal for An) 

N 

'5e 



s^ '^[PcA'^beiow ieft{vi) is X - pivotal for A„) 

i=l 

+ Pe,s{ubelow rightivi) is X ~ pivotal for An)] 
1 - s 1 1 1 " 

< (— — — + — -— + 2s—)y^Pe,sivi is X- pivotal for An). 

0(1 — e) 0(1 — e) 0€ be 



So there is a continuous positive function 7(e, s) such that 

JV N 

Pe,s{vi is u — pivotal for An) < 7(e, s) Pe,s{vi is X — pivotal for An) 

1=1 i=l 

and by using the Russo-like formulas ^ and ^ from above we get ^^(e, s) > 
dj{e, s)^^{e, s). Since j{e, s) is positive and continuous, for any a > we can 
find M and (j) € (0, f) such that M > S^{e,s) on [a, 1 — a]^ and tanc/) — M. 
The directional derivative of 9„(e, s) thus satisfies 

Y7r> ^ ^ ■ A\ ' ■ ■ 

Vfc)„ • (cos 0, — sm (p) = — — cos 



de ds 

dQn , 5e„ . , 

> 07—;; — COS (p — sm ffl 

OS OS 

= — cos 4'~^~ (tan (j) — 6j) 

>0 (3) 

Set a < min{-^ec, ^s) and mark points a, 5, and c (see Fig. [7]) such that a 
and 6 are inside of [a, 1 — a]^, a is at height s but left of Ec, & is such that the 
vector from a in the direction of (cos 0, — sin 0) of appropriately chosen length 
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s _ 



b 




Figure 7: 



crosses the line e — Sc and lands at b, and c has the same e-coordinate as b and 
has s-coordinate equal to 0. From ^ we conclude that 6„(a) < 0„(&) < 6„(c). 
Taking n ^ oo and using that 7(5, s) does not depend on n, 0(e, s):=P^^s{There 
is an infinite path from the origin) satisfies 6(a) < 0(&) < 6(c). Since e > ec 
at c, we have 0(c) =0 and therefore Q{a) = 0. This completes the proof of 
Lemma 2.n 

Proof of Theorem 2: We intend to couple the processes Z„ with two 
different initial color sequences Aq and Bq in the same product space of the 
i.i.d. X{v) and Y{v) random variables to show that for any finite number of 
colors q > 2 and for e in some interval (e^, 1) with e'^ < Cc (so that percolation 
of arrows occurs for e £ (e^,ec)), the Zn process is ergodic. We note that the 
value of obtained from this proof depends on the value of q. 

Let A and B refer to the two processes of Z„ corresponding respectively 
to the initial color sequences Aq and Bq (with A{v) and B{v) referring to the 
random colors of a single vertex v). Let C be the {0, l}^-valued process with 
C{v) = if A{v) = B{v) and C{v) = 1 if A{v) ^ B(v). The goal is to 
show that almost surely C eventually turns to all zeros (in any finite spatial 
interval). To achieve this goal we take advantage of the fact that the state 
space {0, 1}^ is partially ordered and use monotonicity present in the system 
(for more information see |13)). Let vi and Vr be the vertices above left and 
above right of v. We can bound the one step transition probabilities of the C 
process as follows. 



This is so because when C{vi) = C{vr) = 0, then A{vi) = B{vi), A{vr) = B{vr), 
and A{v) must equal B{v). The case C{vi) = and C{vr) = 1 means A{vi) = 



P{C{V) = l\C{vi),C{Vr)) 




1-P(0) 







if C{vi) 
if CM 
if C{vi) 
if Civi) 



C{Vr) = 0, 

0, C{Vr) - 1, 

1, C{Vr) = 0. 
C{Vr) = 1. 
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B{vi) and A{vr) ^ B{vr). To have A{v) ^ B{v) we need "z^" at v or "\/>" 
at w in which case we know that at least one of A{vi) ^ A(vr) or B(v{) ^ B{vr) 
happens (because otherwise we would have C{vi) — C{vr)) and the probability 
that the newly chosen color is different from the color of the other process is 
at most ^ (it is if both A(yi) ^ A{vr) and B{vi) ^ Bivr)). The case 
C{vi) = 1 and C{vr) = is similar. For the case C{vi) = 1 and C{vr) — 1, even 
though it is stih possible that C{v) = when X{v) = "\/^", A{vi) ^ A(ur), 
Bi^vi) 7^ B{vr), and the value of the newly chosen color Y{v) goes to both A{v) 
and to B{v), we bound the transition probability by 1 — P{(D) for simplicity. We 
keep the initial configuration of the C process for a newly defined {0, l}^-valued 
C process to which we assign the transition probabilities 

P{C'{V) = l\C\vihC'{Vr)) 

if C'{vi)^C'ivr)^0, 

P{\)+P{\^)'^ if C'{vi) = l,C'{Vr)^0, 
1 - P(0) if C'ivi) = C'iVr) = 1. 

(4) 

This modification of the transition rule from C to C can be obtained with an 
appropriate coupling in which the only change is that some vertices switch from 
to 1. An appropriate coupling of C and C is the following. The probability 
space for the coupled process will be the product space of the i.i.d. arrow- valued 
random variables X(v) and i.i.d. uniform[0,l] random variabls U{v) at each 
vertex v € V. For the vertices v in the initial row, A{v) and B{v) are assigned 
the values at v of the color sequences Aq and Bq respectively. The {0, l}-valued 
random variables C{v) and C'{v) take the value if A{v) = B{v) and value 1 if 
A{v) ^ B{v). For the subsequent rows, the processes A and B evolve according 
to the rules of the Z„ process with the values for the new random color Y{v) 
chosen according to the value of U{v) by splitting the range of values [0, 1] into q 
equal sub-interval and assigning the q colors to each of the sub-interval in some 
predetermined order. C{v) = if A{v) — B{v) and C{v) = 1 if A{v) ^ B{v) as 
above. For the C process we want to define the transition rules in such a way 
that the statement "if C(w) = 1, then C'{v) = 1" always holds and such that 
the marginal transition probabilities of the C" process are as in 

For V in the initial row, the statement "if C{v) = 1, then C'{v) = 1" holds 
since C{v) — C'{v). At each vertex v after time look at the values of all 
the processes involved at the vertices vi and Vr in the previous row and at the 
values of X{v) and U{v). Suppose that in the row above which contains vi and 
Vr it is true that whenever C{v) — 1, then C'{v) — 1. If X{v) = "0", assign 
C'(v) = C{v) = 0. If X{v) = " \" or X{v) = " /^", then the value for C'{v) is 
chosen to be C'{vi) or C'{vr) respectively, and the statement "if C{v) — 1, then 
C'{v) = 1" carries over from the previous row. In case X{v) = " \/^" and 
C'{vi) and C'{vr) agree in value, C'{v) takes that value. If the values of C'{vi) 
and C'{vr) disagree, the rules for choosing the new value for C'{v) will depend 
on the values of A{vi), A{vr), B{vi), and B{vr). Unless it is a case where a 
special rule is required, choose the new value for C'{v) in the similar way as 
the value for Y{v) is chosen with value assigned to the first sub-interval [0, 
and value 1 assigned to the rest. If C{vi) = C{vr) = 0, then C{v) = and the 
statement holds. If C{vi) = C'{vi) = and C{vr) = C'{vr) = 1, consider the 
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values of A{vi), A{vr), B{vi), and B{vr)- If A{vi) ^ A{vr) and B{vi) ^ B{vr), 
then C{v) = and the statement holds. If A{vi) ^ A{vr) and = B{vr), 

the special rule is required, as follows. Choose a value for Y{v) and a value 
for C'{v) based on the outcome of the uniform[0,l] random variable U{v) the 
following way: split the interval [0, 1] into q equal sub-intervals and assign the q 
colors to them in the same predetermined order as before for the Y{v) value, and 
for the C'{v) value assign to the sub-interval corresponding to the color B{vr) 
and 1 to the other sub-intervals. The special rule is consistent with the rules for 
choosing values at v for the processes A, i?, and C, the transition probabilities 
for C'{v) correspond to ([4]), and the statement "if C{v) — 1, then C'{v) — 1" 
holds. If A{vi) = A{vr) and B{vi) ^ B{vr) the special rule is the same except 
the value for the C'{v) is assigned to the sub- interval that holds the color of 
A{vr). The case when C{vi) — C'{vi) — 1 and C{vr) — C'{vr) = is similar, 
and in all other cases the statement holds trivially. We can say that the C" 
process dominates the C process in the sense that the statement "if C{v) = 1, 
then C'{v) = 1" always holds. 

Therefore, if with any initial configuration the new process C'{v) turns to 
all zeros (or "dies") eventually almost surely (in any finite spatial interval), 
then in the original color and arrow process Zn the colors of all the vertices (in 
any finite spatial interval) will eventually become the same for all initial color 
configurations, which implies the uniqueness of the invariant distribution and 
convergence to it. 

To show that the C process eventually turns to all zeros, we use the auxiliary 
enhancement result of Lemma 2. To relate the random graph G of X(v) arrows 
to a {0, 1}^ valued process, we think of all the vertices of G that have an infinite 
path coming out of them as having value 1 and all the vertices whose (directed) 
cluster of arrows is finite as having value 0. With this in mind, we define a new 
{0, 1}^ valued process C" by assigning C"{v) = 1 to all vertices at the initial 
time, by directing time down (opposite of arrow direction), and by transfering 
value of 1 down by any arrow pointing to a site with value 1. The transition 
probabilities of C" are the following: 



An almost sure limiting configuration of all zeros of the C" process is equivalent 
to non-percolation of arrows. The enhancement function described in the arrow 
percolation context is activated with probability s at a vertex v only when we 
have at V a double-arrow and at least one of the vertices above right or above left 
of V has "0" value. For the C" realization with an X{v) arrow environment to 
activate the enhancement it is necessary that at least one of the vertices above 
right or above left has value 0. But that is not sufficient since it is possible that 
both vi and Vr have arrows coming out of them but with all the paths coming out 
of, for example, vi terminating at "0" points. If we denote the {0, l}^-valued 
enhanced percolation of arrows process by C*, its transition probabilities will 



P{C"{V) = l\C"{vi)fi"{Vr)) 






P{\) + P{\/^) 
1 - P(0) 



if C"{vi) 
C"{vi) 
C"{vi) 
C"{vi) 



C"(«,) -0, 

0, C"(t;,) = 1, 

1, C"{Vr)=0, 
C"{Vr) = 1. 
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satisfy the following: 



P{C*{v)^l\C*{vi),C*{Vr)) 

if C*ivi)=C*ivr)^0, 

P{^) + P{\/'){1 - S) if C*{vi) = 0, C*{Vr) - 1, 

^(\) + - S) if C*{VI) = 1, C*{Vr) = 0, 

1-P(0) tf C*{vi)^C*{Vr) = l. 



> 



Set s — - > 0. For that choice of s with the standard coupling we can dominate 
the process C by the process C* in the sense that "if C'{v) = 1, then C*(v) = 
1" . Lemma 2 says that e^"''=e^ < e^- This means that for the interval of 
e'^ < e < 1 the enhanced cluster of the origin is almost surely finite, and all the 
values come from "0" points which give value 0. Therefore, the unique invariant 
distribution for such e for the processes C*, C , and C is the delta measure at 
all zeros. □ 



Remark: The dominating C process considered in the proof of Theorem 2 
is a simpler {0, l}^-valued process for which many techniques are developed in 
discrete and continuous time ([71[Tni[I3] among others). For example. Theorem 
4.1 of Chapter I or the proof of Theorem 3.32 of Chapter III of [TB] provide for 
an easily derived lower bound such that for e > the continuous time version 
of the C" process a.s. converges to the all zero configuration from any initial 
state. The proof of Theorem 3.32 easily extends to the discrete time version of 
C providing the same lower bound. However, we chose to use the enhancement 
argument of Lemma 2 because it allows us to compare the actual critical values 
tc for percolation and for ergodicity. 



Proof of Theorem 3: Let V C V he the set of vertices from which only 
the horizontal half-line to the right of and including the origin can be reached by 
arrows taken from the set of all possible arrow configurations {\, \/^, 0} 
at each vertex v £ V. See Fig[S] 

aaaaaaaaaa 



• 



Figure 8: Set V. 

On the space = {\,/^,\/^, 0}^ x {a, &, c, ...}^, where a,b,... are q 
colors, let P be the product distribution of the i.i.d. arrow- valued X{v) and 
the color-valued uniform Y{v) random variables for v E V. For each lj E 
and for each v E V define new arrow- valued random variables W{v){u!) and 
color- valued random variables Z{v)(ll)) as follows: for v in the top row of V' fix 
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W{v){uj) — "0" and assign some configuration of colors Z{v) to the top row. 
For tlie vertices in each, consecutive row below the top one define 



W{v){uj) = 



\ if X{v){lo) = \, 

^ if X{v){lo) = /, 

\^ if X{v){u:) andZ{vi){uj)^Z{vr){uj), 

otherwise, 



(here vi means the vertex above left of v and Vr the vertex above right of v) and 
Z{v){u) = 



Z{vi){uj) if W{v){lo) = \, 

Z{Vr){uj) if W(V){U) - /, 

Z{vi){u) = Z{vr){u^) if W(v){u:) 

Y(v){u) if W(v){u) = 0. 



Thus defined, Z{v) is the same as the color of the vertex v of the Z„ process. 
We want to estimate how the probability that the color Z{v) comes from the 
top row changes. Let G be the collection of vertices v £V' that are connected 
by paths of PF-arrows to the top row. Let vq,vi,... be the leftmost vertices of 
V as in Fig. [8] and consider the sum 

oo 

4=0 

^^ r., ^,P{vieG) ^, ^,P{vi gG) P{v2eG) 
Pivo e G) + P(.o e G)j^^^ + P(.o e + .... 

Denote the upper right vertex of Wi+i by Ui and the event {At least one of Vi 
or Ui is in G and the other one has the same color} by Ei. Then 

P(w,+i e G) =P(X(w,+i) =\ )P{v, e G) + P{X{v,+i) )P(u, e G) 

+ P(X(z;,+i)=\^)P(£;,). (5) 

First we consider the case where all the colors in the top row are color a. 
Then 

P{Ei) ^P{vi e G and Ui e G) + P(wj G G and Ui ^ G and Z{ui) = a) 
+ P{ui e G and Vi ^ G and Z(vi) = a). 

If Ui ^ G, the color of Ui was determined at some vertex v with W{v)="Q" by 
choosing a color uniformly out of the q possible ones. Therefore 

P{vi e G and Ui ^ G and Z{ui) = a) = —P{vi G G and Ui ^ G) 

= P(ui € G and Vi ^ G and Z{vi) — a), 

and so 

P{Ei) = P(i;, e G and Ui G G) + |p(w, G G and Ui ^ G). 

When q> 2, we thus have P{Ei) < P{vi G G) and so 
Piv,+i G G) 

< [P{X{v,+i) -\ ) + P(XK+i) ) + P(X(«,+i) )]P{tH G G) 

= (1 - <5e)P(w, G G). 
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Thus 

P{v^+l € G) 

P{v^ e G) 



<1-Se<l \/i 



when S > and e > 0. Therefore, X^i^o Pi'^i G G) < oo. By the Borel-Cantelh 
lemma P(G contains infinitely many Vi) — 0, or almost surely 3N such that 
Vi >N,v^^ G. 

For any other than the constant a configuration of colors at the top, we can 
compute similarly 

P{Ei) = P{vi £ G and Ui € G and Z[ui) = Z{vi)) 
+ P{vi e G and Ui ^ G and Z(ui) — Z{vi)) 
+ P{ui G G and Vi ^ G and Z(vi) — Z{ui)) 

< P{vi e G and Ui e G) + P{vi e G and Ui ^ G)- + P{ui e G and v.^ ^ G)- 

< P{vi e G and e G) + P{vi e G and <^ G)i + P(?ij G G and v., ^ G)i 

= P{v, e G)\ + P{u^ e G)i < sup Piv e G). 

^ ^ V in row i 

We next use ^ and the reasoning in the paragraph after ^ to conclude first 
that for any vertex in row i + 1, P{vi e G), P{ui S G), and P{Ei) are all 
not greater than 

sup P{v e G), 

i; m row i 

and second that for all the vertices u^+i in row i + 1 

P{v.,+i e G) < (1 - Se) sup P(w e G). 



V in row t 



Hence, 



and 



sup P{v e G) < (1 - 6e) sup P{v e G), 

t; in roio i-^-l v in row i 



J2 ^(«« e G) < ^ sup P{v e G) < cx). 

Therefore, for any initial coloring, almost surely Vi is not connected by the W- 
arrows to the top of V for large enough i. The same is true for the W^-arrow 
cluster of any deterministic finite interval of the i*'* row of vertices. We can 
conclude that for any subsequence limit distribution, any joint distribution of 
colors of a finite collection of vertices is color permutation invariant. □ 



Proof of Theorem 4: We couple the q = oo process with any initial 
color configuration (or equivalence class partition) to the one with the initial 
row of a constant color. By the argument of the proof of Theorem 3, for any 
finite interval of vertices, if we start with the constant color configuration some 
sufficiently large amount of time earlier, the cluster G of vertices that can be 
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connected by the W{i;) arrows to the initial row will not intersect the finite 
interval. For g = oo, we show that colors of the vertices not in the cluster 
G obtained with the constant color initial configuration are almost surely the 
same for all initial color configurations. This will follow if we can show that 
W{v) = "0" with the initial constant color implies W{v) = "0" for any other 
initial color configuration because at such vertices a new color is chosen which 
is independent of the initial configuration. And this is true because, given that 
the initial colors are constant, W{v) = "0" happens when either X{v) = "0" 
and, therefore, W{v) = "0" for any initial color configuration, or when X{v) = 
"\/"' and the colors of vi and Vr arc different. If we start with the constant 
color configuration, for vi and Vr to be of different color at least one of them, say 
Vr, has to be connected by the X-arrows to some vertex w that is located later 
than the initial time and such that X{w) ~ "0" (or both, vi and u^, have to 
be connected to distinct such w). For the case q = oo (when every "0" vertex 
almost surely creates a new color that does not match any previously existing 
color), every '"^y^" along the path of the X{v) arrows from Vr to w will be 
changed to "0" (unless the double arrow creates branches that come together 
again before reaching w). But in this situation the same thing will happen 
along the path from v to w for any initial color configuration, and colors of the 
vertices above left and above right of v will also be different no matter what 
the initial color configuration was. Therefore W{v) = "0" for every initial color 
configuration. Therefore, if starting from the constant color initial configuration 
the distribution of equivalence class partition at time n converges as n — oo, 
then the q = oo system is ergodic. 

The fact that the limit distribution of the equivalence class partition exists 
as n — >■ oo for a constant color initial configuration follows from the argument 
above because it allows us to couple equivalence class partitions corresponding 
to different initial times. Instead of looking at the distributions of the equiv- 
alence class partitions of a finite interval of vertices at time n and time n + 1 
corresponding to the constant color initial configuration, we look at the interval 
at time and compare the distributions of its equivalence class partitions corre- 
sponding to the initial times —n and —n—1. The equivalence class partition of 
the interval corresponding to the initial time —n—1 can be coupled to the one 
corresponding to the initial time — n by adding a row of the X arrows at time 
—n that point to time —n— 1 and moving the initial constant color configuration 
from row —n to row —n—1. By the argument above, if n is large enough, colors 
of the interval with high probability will be the same for all initial configura- 
tions at time —n as they would be for the constant color initial configuration. 
Therefore, moving the constant color initial configuration from time —n to time 
—n — 1 will change colors at time —n but not the colors of the interval. □ 
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